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Abstract 

It is an elementary fact that the size of an orthogonal array of strength t on k factors 
must be a multiple of a certain number, say Lf . that depends on the orders of the factors. 
Thus Lt is a lower bound on the size of arrays of strength t on those factors, and is no 
larger than L the size of the complete factorial design. We investigate the relationship 
between the numbers Lf, and two questions in particular: For what t is Lt < L *.? And 
when Lt = L is the complete factorial design the only array of that size and strength 
t? Arrays are assumed to be mixed-level. 

We refer to an array of size less than L & as a proper fraction. Guided by our main 
result, we construct a variety of mixed-level proper fractions of strength k — 1 that also 
satisfy a certain group-theoretic condition. 
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1 Introduction 


Let D be an orthogonal array of size N on k factors, the ith factor having s t values or levels. 
It is easy to see that if t is the strength of D , then N is a multiple of 


L t 




( 1 ) 


where |/| is the number of indices in /. These points are reviewed in detail below. 

Obviously L t is a lower bound for N , although an array of strength t and size N — L t 
may not exist. Moreover, L t < Si---Sk = Lk, the size of the complete factorial design on 
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these factors, since the complete design certainly has strength t. From the point of view of 
applications, it is of interest to know when L t < L k , since in this case one may seek a design 
that is a proper fraction of the full factorial design. Theorem 12.11 gives a simple criterion based 
on the orders s* to determine those t for which we have L t < L k . 

1.1 Definitions and notation 

As indicated above, we use |/| to denote the cardinality of the set I. 

A complete factorial design is a Cartesian product A\ x • • • x A k , where A; is a finite set. 
In statistical design, Ai is the set of levels of factor i , and the elements of A\ x • • ■ x Ak are 
treatment combinations. We will refer to s; = \Af\ as the order of factor i. An orthogonal array 
or design D on these factors is a multisubset (a subset with possible repetitions) of A\ x • • • x A k . 
The size of the array, N, is the number of elements (or runs), counting multiplicities. The 
design is symmetric if si = • • • = Sk, and otherwise is asymmetric or mixed-level. We will say 
that a design is a proper fraction if it is a a proper subset of a complete factorial design. 

If we write the elements of a design D as columns, then we may represent D as a k x N 
matrix of symbols. The projection of D on j factors is then the j x N submatrix consisting 
of the rows corresponding to those factors. We say that D has strength t if for any subset 
/ C {1,..., k} of size t, the projection of D on the factors indexed by / consists of A j copies 
of the Cartesian product riie/ Ai, for some A j. Evidently we have 

n=\j 

iei 

from which it follows easily that N is a multiple of the number L t given by (JTJ) . In particular, 
Li = lcm(s'i ,... ,s k ) and L k = s i • • • s k , and it is easy to see that 

Lt\Lt+i 

for each t — 1,..., k — 1 (each product in L t is contained in a product in L t+1 ). Therefore 

L\ < L 2 < • ■ ■ < L k . (2) 

In a symmetric design with si = • • • = s k = s we have L t = s t . 

For each t, L t is the smallest possible size of an array of strength t. An array of this size 
may not actually exist. For example, in a symmetric design on k — 4 factors, each with 6 
levels, we have L 2 = 36, but there is no array of strength 2 and size 36, since this would be 
equivalent to two mutually orthogonal Latin squares of order 6 (see, e.g., [3], pages 11 and 33). 

As indicated above, it is of interest to know for which t we have L t < L k . Theorem 12.11 
gives a useful criterion, and also strengths the inequality ((2]). 

2 The numbers L t 

For each / C {1,..., k} let 

ei = gcd (s i} i £ I), 
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and consider those / for which ej > 1. Let 


d = max{|/| : ej > 1}. 

Thus |/| > d implies ej = 1. We see that 1 < d < k, and d = 1 iff the orders s, are pairwise 
relatively prime. At the other extreme, d = k iff gcd(si,... ,Sk) > 1. For example, d — k if 
the design is symmetric. 

In the proof below we need to choose a subset / of size d for which ej > 1. There can be 
more than one such subset: for example, if the numbers are 8,12,18, 27, then d — 3, since 
the four numbers have no common factor but {8,12,18} and (12,18, 27} have gcd greater than 
1; the corresponding sets / C (1,..., 4} are (1, 2, 3} and {2, 3,4}. 

Theorem 2.1. We have L\ < ■ ■ ■ < L c i = • • • = Lk- In particular, L t = Lk iff t > d. 

Before proving this theorem, we mention two useful facts. First, if Si,, S n are sets of 
positive integers and S = U*5}, then 


lcm (S') = lcm(lcm(Si),..., lcm(S n )). 


(3) 


Second, for positive integers ai,..., a n we have 



^1 ’ * * ^n 

gcd(ai, ...,a n )' 


(4) 


These can both be proved prime-by-prime: Define ord p (6), the order of the prime p in b, to be 
the largest power of p dividing b. To prove (J4j), for example, let fi = ord p (aj). Then (J4J) is the 
statement that for each prime p, 


( n \ n 

Efi-fj min fr 

i =1 / i =1 J 

(Property flD) is familiar in the case n — 2, namely, ab = lcm(a, b) gcd(a, b).) 

Proof of Theorem \2.1[ Fix t. In order to compare L t and L t+ 1 , we begin by organizing the 
subsets / C {l,...,fc} of size t into overlapping families Bj indexed by the subsets J of size 
t + 1. Namely, we put 


Bj = {I C J: \I\ = t}. 
Let B = {I C {1,..., k }: |/| = t}. Then 

b = U B -” 
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Now from (J3J) and (JU) we have 



( 5 ) 


If t > d then the expression ([5]) is exactly L t +±. For if t > d then ej = 1 whenever 
| J\ — t + 1. Tims t > d implies that L t = L t+ \. 

If t < d, we claim that (J5]) is strictly less than L t+ \. To see this, fix a prime p that divides 
exactly d of the numbers si,..., s (This is possible as there exists a set / C {1,..., k} with 
\I\ — d and ej > 1; let p|e/.) Our claim follows as soon as we show that the order of p in ((5]) 
is strictly less than ord p (L t+ i). To do this, we need to show that this holds for each term in 
d5]), that is, 

ordp (n Si/ejj < ord p (L m ) (6) 

for each Jc{l,...,fc} with | J\ — t + 1. 

Let fi = ord p (si). Note that exactly d of the integers /* are positive. Renumbering if 
necessary, we may assume without loss of generality that 


fi > /2 > • • • > fd. > 0 — fd+l — ■ ■ ■ — fk- 


Now on the one hand, the order of p in L t +i is the maximum of the orders of p in the 
products riiej \J\ = t + 1. But this maximum is attained for J — {1,..., t + 1}, and equals 
fi H-1- ft+i, that is, 

t +i 

ord p (L i+1 ) = (7) 

i —1 

On the other hand, for any J C {1,..., k} we have 

ordp(ej) = min /•, 
ieJ 

and so 

or dp JJ Si/ej = /< - min f t . 

\i( J J H J 

We focus on sets J with | J\ — t + 1, and consider two cases: 
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• mirijgj f, > 0: Here 


t +1 

Z fl ~ fi < Z fi ~ Z U 

i£j zEJ i= 1 

• mirijgj /,; = 0: In this case, f % — 0 for at least one i 6 J, so 

_ _ t+i 

Z ^ _ ™ ■/' = Z fi < Z fi - 

i£j i&J i =1 

In either case (EJ) holds, and therefore the order of p in (J5J) is less than (JTJ), which is what we 

needed to prove. □ 

Theorem 12.11 provides useful information about possible orthogonal arrays without our 
directly having to calculate the numbers L t . In particular, it gives a necessary condition for 
the existence of a proper fraction of given strength, as illustrated by the examples below. A 
special case of the theorem is given in [2], 

Example 1. In a 2*3 J experiment, that is, one having i + j factors of which i have 2 levels and 
j have 3, we see that d = ma x(i,j). In this case, no proper fraction has strength t > ma 

Example 2. Consider a2x3x5x6xl0xl5 experiment. Any subset of {2, 3, 5, 6,10,15} of 
size 4 must have ej = 1, as two of its elements must be relatively prime. On the other hand, 
there are subsets of size 3 that have ej > 1 - for example, {2,6,10} (that is, / = {1,4,5}). 
Thus d — 3, and so no proper fraction has strength 3. 

Example 3. In order for there to exist a proper fraction of strength k — 1 in an experiment 
with k factors, it is necessary that d = k. That is, the orders S\,... ,Sk must share a common 
factor. In Section [4] we give methods of construction that generate a number of proper fractions 
of strength k — 1 that also satisfy a group-theoretic condition when the set of treatment 
combinations is a nonabelian group. The condition and its origin are discussed there. 


3 Uniqueness 

When L t = L^, the smallest array of strength t has the same size as the complete factorial 
design, and it is natural to ask whether the complete factorial is the only array of its size of 
strength t. It is easy to see that for strength t — 1 this is not true. For suppose factor i is 
indexed by the set A* where |A,| = Sj. Fill a kx matrix by putting Lk/si copies of Ai in row 
i, in an arbitrary order. The resulting array has strength t — 1 and size Lk, but it is easy to 
fill the rows in such a way that the columns do not consist of all the elements of A\ x • ■ • x Ak- 
Of course, the rows of an array of strength 1 can be filled independently of each other, 
and so one might conjecture that such a construction is impossible if we require strength 2 or 
higher, but in fact this is not true either. Consider a 3 x 2 x 2 factorial experiment and let 
A\ = {0,1,2} and A 2 = A 3 = {0,1}. In this case L 2 = 12, the size of the complete factorial 
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design, and so there is no smaller orthogonal array of size 12 and strength 2 on these factors. 
We easily check that the 3 x 12 array 

/ 001122001122 \ 

000000111111 

yoooiiiiioioo/ 

has strength 2, but it is clearly not the complete design, as some ordered triples occur more than 
once while others are missing. It would be interesting to know whether there are applications 
in which such an array might be a useful alternative to the full factorial design. 

4 Constructions 

In this section we construct a number of mixed-level orthogonal arrays of strength k — 1 on 
k factors. They are all proper fractions and in addition satisfy a group-theoretic property 
(a “conjugacy” condition) given in [TJ that we describe below. As we saw in Example [3j 
Theorem 12. II requires that the orders s* share a common factor. 

Our examples consist of k = 3 to 6 factors, with orders s t < 10. Verification that each array 
satisfies our two conditions (strength k — 1 and conjugacy) was done by computer, for which 
the code can be found in [2]. Only those cases listed in Table [1] (below) have been checked. 

4.1 Conjugacy 

Throughout this section we will suppose that the sets Ai indexing the levels of the factors are 
groups, and will use Gi rather than A, as a more suggestive notation. The set G = G\ x • • • x G k 
of treatment combinations is a group, and is abelian iff all the factors Gi are abelian. 

Recall that the elements a and b of G are conjugate if b = gag for some g e G. Conjugacy 
is an equivalence relation on G, and so partitions G into conjugacy classes. In an abelian group, 
the conjugacy classes are all singleton sets. It is an elementary fact that the conjugacy classes 
of G\ x • • • x Gk are of the form Cj x ■ • ■ x C k where Gi is a conjugacy class of Gi. The condition 
we require is that the design D be a union of conjugacy classes ofG$ 

In all our examples G i, and therefore G, is a nonabelian group, while Gi will be abelian 
for i > 1. Thus the conjugacy classes of G are essentially those of G\ - namely, of the form 
C x {g 2 } x • ■ • x {g k } where C is a conjugacy class of G i and g % G Gi. The nonabelian groups 
we will use are the following, where e will denote the identity element. 

• S 3 , the symmetric group on 3 letters. We have S 3 = {e,x,y,a,b,c}, where x and y are 
the 3-cycles and a, b and c are the 2-cycles. We will use either of two orderings: 

1. e\xy\abc. 

lr The condition assumed in [Tj is that the counting function of the multiset D be constant on conjugacy 
classes. Since we are constructing proper fractions, D is a set and its counting function is its ordinary indicator 
function, which is constant on conjugacy classes iff D is a union of such classes. 
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2 . e | a b c \ x y. 

The vertical lines are inserted merely to indicate the conjugacy classes. 

• Dih 4 = {e, q, r, s, a, b, x, y}, the symmetries of the square, where q is the half-turn, r and 
s are the quarter-turns, a and b are the reflections about the diagonals, and x and y are 
the reflections about the lines joining the midpoints of the opposite sides. The conjugacy 
classes are e\q\r s \ a b\xy. 

• Dih 5 = {e, a, b , c, d, v , w, x, y, z }, the symmetries of the regular pentagon, where a and d 
are the rotations of ±72°, b and c are the rotations of ±144°, and v, w, x, y and z are the 
reflections about a line from a vertex to the midpoint of the opposite side. The conjugacy 
classes are e\ad\bc\vwxyz. 

In each of the above, the ordering of conjugacy classes is fixed, and the ordering of elements 
within each class is fixed but arbitrary. These fixed orderings are assumed in the constructions 
below. In general, the group Dih n is a dihedral group, and represents the symmetries of the 
regular n- gon. Of course, S 3 is also Dihz. 

For abelian groups we will use Z n , the integers modulo n, whose elements are 0,1,..., n. 
We will fix this order. 

4.2 Construction 

If Si = 6 , 8 or 10, we let G 1 = S3, Dih 4 or Dih 5 , respectively. For i > 2, we let G t = 7L Si . Our 
construction has the following cases. 

gcd(si,..., Sfc) = 2 or 4: We construct an array of minimal size Defining the 

integers 

Lk -1 
v\ = -, 

Sl 

Vj = ^ k '~ 1 for j > 2 , 

s 2 • • • Sj 

we fill in the rows of the array as follows: 

Row 1: Write the elements of G\ in the fixed order, the whole sequence repeated V\ times. 
If Si = 6 , use S3 in its first ordering. 

Row 2 : Write the elements of G 2 (= Z S2 ) in the fixed order, repeating each element u 2 times. 
Rows 3 through k — 1: Write the elements of Gj[= Z s .) in the fixed order, repeating each 
element Vj times. Repeat the whole pattern until the row is complete (a total of L^_ 1 entries). 

Row k: Write the elements of Gk, repeating each element Vk times, in the following pattern: 
first in the given order, then in the reverse order, alternating in that way until the row is 
complete. 

Example 4. The following arrays illustrate the construction method we have described. 
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1/2 fraction ofa6x2x2x2 design, strength 3. Here G = S 3 x Z 2 x Z 2 x Z 2 . 


/exy a bcexy a bcexy a bcexy a bc\ 

000000000000111111111111 
000000111111000000111111 
\ 0 0 0 1 1 1 1 1 1 0 0 0 0 0 0 1 1 1 1 1 1 0 0 0 / 


• A 1/2 fraction of an 8 x 2 x 2 design, strength 2 . G — Dih± x Z 2 x Z 2 . 


( e q r s a b x y e q r s a b x y\ 

0000000011111111 

0000111111110000 / 

• A 1/2 fraction of a 10 x 2 x 2 design, strength 2 . G = Dih 5 x Z 2 x Z 2 


( e a b c d v w x y z e a b c d v w x y z\ 

00000000001111111111 
0000011 1111111100000 / 

• A 1/4 fraction of an 8 x 4 x 4 design, strength 2 . G — Dih 4 x Z 4 x Z 4 . Vertical lines are 
inserted to reveal the pattern in the last row. 


( e 

Q 

r 

s 

a 

b 

X 

y 

e 

Q 

r 

.s 

a 

b 

X 

y 

e 

Q 

r 

s 

a 

b 

X 

y 

e 

9 

r 

s 

a 

b 


y\ 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

3 

3 

\0 

0 

1 

1 

2 

2 

3 

3 

3 

3 

2 

2 

1 

1 

0 

0 

0 

0 

1 

1 

2 

2 

3 

3 

3 

3 

2 

2 

1 

1 

0 

0 / 


gcd(si,..., s*;) = 6: We in fact consider only symmetric “6 fc ” designs, so that Lk-i = 6 fe_1 . 
The method described above may be used to create arrays of this minimum size (and therefore 
1/6 fractions), but the arrays will not satisfy the conjugacy requirement. If we modify the 
integers as follows: 


ni = —— = 3 • 6 fc_2 , 

Vj = —— k 1 = 3 • for j > 2, 

s 2 • • • Sj 

then the method will produce 1/2 fractions of strength k — 1. Note that the first ry = 3 
elements of G'i = A 3 is a union of conjugacy classes of S 3 . 

We don’t present an example since N = (l/2)6 fc rather large. 

gcd(si,..., Sfc) = 3: We apply this to 6 x 3 x 3 • • • factorial experiments. Here we alter 
both the integers v t and the steps of construction, since the original steps will produce arrays 



of minimum size (1/3 fractions) that do not satisfy the conjugacy property. We set 

2-hfc-l _ 0 k-2 

Vi = -= 2-3 , 

Sl 

= 2 ^ : 1 = 4 .for j > 2 , 

s 2 • • • Sj 

and fill in the rows of the array as follows: 

Row 1: We use the second order of S 3 . Write the elements of G\ = S 3 in this fixed order, 
then in the reverse order, alternating until there are N = 2Lk-i- 

Rows 2 through k — 1: These steps are identical to those given above. 

Row k : Write the elements of Gk = Z 3 , repeating each element Vk = 4 times. Then do the 
same, but permuting the elements of Z 3 cyclically, and again with another cyclic permutation, 
continuing in this pattern until the row is filled. 

This method produces 2/3 fractions. 

Example 5. A 2/3 fraction of a 6 x 3 x 3 design, strength 2. Here G = S 3 x Z 3 x Z 3 . Vertical 
lines are inserted to reveal the pattern in the last row. 


( e 

a 

b 

c 

X 

y 

y 

X 

c 

b 

a 

e 

e 

a 

b 

c 

X 

y 

y 

X 

c 

b 

a 

e 

e 

a 

b 

c 

X 

y 

y 

X 

c 

b 

a 

e \ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

\0 

0 

0 

0 

1 

1 

1 

1 

2 

2 

2 

2 

2 

2 

2 

2 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

0 

0 

0 

0 / 


Table |T] summarizes 31 arrays constructed using the methods we have described. 
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Table 1: This table lists arrays constructed using the methods in this paper. Strength is A; — 1 
where k is the number of factors. All arrays are if minimal size L k _i unless otherwise indicated. 
Examples are separated into four groups according to gcd(si, ..., s k ). 


Complete 

Size of 

Size of 



Design 

Complete Design 

Array 


Fraction 

6x2x2* 

24 

12 


1/2 

6 x 2 x 2 x 2* 

48 

24 


1/2 

6x4x4 

96 

48 


1/2 

6 x 4 x 4 x 4 

384 

192 


1/2 

6x4x2 

48 

24 


1/2 

6x6x2 

72 

36 


1/2 

6x6x4 

144 

72 


1/2 

8x2x2 

32 

16 


1/2 

8 x 2 x 2 x 2 

64 

32 


1/2 

8x2x2x2x2 

128 

64 


1/2 

8x2x2x2x2x2 

256 

128 


1/2 

8x6x6 

288 

144 


1/2 

8 x 6 x 6 x 6 

1,728 

864 


1/2 

8x4x2 

64 

32 


1/2 

8x6x2 

96 

48 


1/2 

8x6x4 

192 

96 


1/2 

10 x 2 x 2 

40 

20 


1/2 

10 x 2 x 2 x 2 

80 

40 


1/2 

10 x 4 x 4 

160 

80 


1/2 

10 x 4 x 4 x 4 

640 

320 


1/2 

10 x 6 x 6 

360 

180 


1/2 

10 x 6 x 6 x 6 

2,160 

1,080 


1/2 

10 x 4 x 2 

80 

40 


1/2 

10 x 6 x 2 

120 

60 


1/2 

10 x 6 x 4 

240 

120 


1/2 

8x4x4 

128 

32 


1/4 

8 x 4 x 4 x 4 

512 

128 


1/4 

6x6x6 

216 

108 

= 3L k _i 

1/2 

6 x 6 x 6 x 6 

1,296 

648 

= 3 L k i 

1/2 

6x3x3 

54 

36 

= 2 L k -i 

2/3 

6 x 3 x 3 x 3 

162 

108 

= 2 Lk-i 

2/3 


* Another example is constructed in [lj. 
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